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Abstract 

We review the current status of theoretical calculations of the hadronic light- by- light scattering 
contribution to the muon anomalous magnetic moment. Different approaches and related issues 
such as OPE constraints and large breaking of chiral symmetry are discussed. Combining results 
of different models with educated guesses on the errors we come to the estimate 

a HLbL = ( 10 .5±2.6) x 1CT 10 . 



The text is prepared as a contribution to the Glasgow White Paper on the present status of the 
Muon Anomalous Magnetic Moment. 



1. Introduction. 

From a theoretical point of view the hadronic light -by-light scattering (HLbL) contribution to the 
muon magnetic moment is described by the vertex function (see Fig. 1 below): 

r (fl), x ■ 6 f rf4fc l f rf4fc 2 ng^g, fci,fc 3 ,fc 2 ) v .. ■ y p , a . . 

r), j (p 2 ,Pi) = te Jj^Jj^y: Y(A-h~m^) 7,(1) 

where is the muon mass and Hl a , pcr (q : ki, k 3 , k^), with g = p 2 — p\ = — k\ — k2 — k 3 , denotes the 
off-shell photon-photon scattering amplitude induced by hadrons, 

n ^p C r('?: fc i> fc 3,fe) = Jd 4 x 1 Jd i x 2 Jd i x 3 exp[-i(fci • xy + k 2 ■ x 2 + k 3 ■ x 3 )] 

x (0|T{j M (0) j„(xi) Jp (x 2 ) j„(x 3 )}\0) ■ (2) 

Here is the Standard Model electromagnetic current, j^x) = Q ' q q(x)"f ^qi^x) , where Q q denotes 
the electric charge of quark q. The external photon with momentum q represents the magnetic field. 




Figure 1: Hadronic light-by-light scattering contribution. 



We are interested in the limit q — > where the current conservation implies that T M 
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is linear in q, 
(3) 



The muon anomaly can then be extracted as follows 
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Unlike the case of the hadronic vacuum polarization (HVP) contribution, there is no direct ex- 
perimental input for the hadronic light-by-light scattering (HLbL) so one has to rely on theoretical 
approaches. Let us start with the massive quark loop contribution which is known analytically, 



HLbL 



(quark loop) = I — I N C Q 
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where N c is the number of colors and m q 3> is implied. It gives a reliable result for the heavy 
quarks c,b,t with m q 3> Aqcd- Numerically, however, heavy quarks do not contribute much. For 
the c quark, with m c ~ 1.5 GeV, 



HLbL 



(c) =0.23x 1CT 10 . 



(6) 
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To get a very rough estimate for the light quarks u, d, s let us use a constituent mass of 300 MeV for 
m q . This gives a HLbL (u, d, s) — 6.4 x 10~ 10 . QCD tells us that the quark loop should be accurate 
in describing large virtual momenta, fcj 3> Aqcd, i-e. short-distances. What is certainly missing 
in this constituent quark loop estimate, however, is the low-momenta piece dominated by a neutral 
pion-exchange in the light -by-light scattering. Adding up this contribution, discussed in more detail 
below, approximately doubles the estimate to ct HLbL w 12 x 10~ 10 . While the ballpark of the effect is 
given by this rough estimate, a more refined analysis is needed to get its magnitude and evaluate the 
accuracy. Details and comparison of different contributions will be discussed below, but it is already 
interesting to point out that all existing calculations fall into a range: 

a HLbL = (11 ±4) x 10~ 10 , (7) 

compatible with this rough estimate. The dispersion of the a HLbL results in the literature is not too 
bad when compared with the present experimental accuracy of 6.3 x 10 -10 . However the proposed new 
.9^,-2 experiment sets a goal of 1.4 x 10~ 10 for the error, which calls for a considerable improvement 
in the theoretical calculations as well. We believe that theory is up to this challenge; a further use of 
theoretical and experimental constraints could result in reaching such accuracy soon enough. 

The history of the evaluation of the hadronic light-by-light scattering contribution is a long one 
which can be found in the successive review articles on the subject. In fact, but for the sign error in 
the neutral pion exchange discovered in 2002 [H [2] , the theoretical predictions for a HLbL have been 
relatively stable over more than ten years. 

Here we are interested in highlighting the generic properties of QCD relevant to the evaluation of 
Eq. Q , as well as their connection with the most recent model dependent estimates which have been 
made so far. 



2. QCD in the Large N c and Chiral Limits 

For the light quark components in the electromagnetic current (q — u ,d , s) the integration of the 
light-by-light scattering over virtual momenta in Eq. is convergent at characteristic hadronic 
scales. We choose the mass of the p meson m p to represent that scale. Of course, hadronic physics 
at such momenta is non-perturbative and the first question to address is what theoretical parameters 
can be used to define an expansion. Two possibilities are: the large number of colors, 1/N C <C 1, and 
the smallness of the chiral symmetry breaking, mj/m^ <C 1. Their relevance can be seen from the 
expansion of a HLbL as a power series in these parameters, 
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where > m M is implied. Only the power dependencies are shown; possible chiral logarithms, 
h^mp/m^), are included into the coefficients c,-. 



Terms leading in the large iV c limit 

The first term, linear in N c , comes from the one-particle exchange of a meson M in the HLbL 
amplitude, see Fig. 2(a). In principle, the meson M is any neutral, C-even meson. In particular this 
includes pseudoscalar mesons 7r°, 77, r/'; scalars /o, &o; vectors n®; pseudovectors a®, /1, /*; spin 2 
tensor and pseudotensor mesons $2, rj2, 1^2 ■ 

The neutral pion exchange is special because of the Goldstone nature of the pion; its mass is much 
smaller than the hadronic scale m p . In a HLbL (7r°) this leads to an additional enhancement by two 
powers of a chiral logarithm [5] , 

Here the 7r°77 coupling is fixed by the Adler-Bell-Jackiw anomaly in terms of the pion decay constant 
w 92 MeV. This constant is O (a/Nc) , therefore N c /F% behaves as a constant in the large- N c 
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(a) (b) 

Figure 2: Diagrams for HLbL: (a) meson exchanges, (b) the charged pion loop, the blob denotes the full 
7*7* — » 7r + 7r~ amplitude. 

limit . The mass of the p plays the role of an ultraviolet scale in the integration over ki in Eq. (j4]) 
while the pion mass provides the infrared scale. Of course, the muon mass is also important at low 
momenta but one can keep the ratio m^/rn^ fixed in the chiral limit. 

Equation @ provides the result for a HLbL for the term leading in the 1/N C expansion in the 
chiral limit where the pion mass is much less than the next hadronic scale. In this limit the dominant 
neutral pion exchange produces the characteristic universal double logarithmic behavior with the exact 
coefficient given in Eq. (0). Testing this limit was particularly useful in fixing the sign of the neutral 
pion exchange. 

Although the coefficient of the In 2 (m p /m ff ) term in Eq. © is unambiguous, the coefficient of the 
]n(mp/m T ) term depends on low-energy constants which are difficult to extract from experiment [2 [3] 
(they require a detailed knowledge of the ir° — > e+e~ decay rate with inclusion of radiative corrections). 
Model dependent estimates of the single logarithmic term as well as the constant term show that these 
terms are not suppressed. It means that we cannot rely on chiral perturbation theory and have to 
adopt a dynamical framework which takes into account explicitly the heavier meson exchanges as well. 

Note that the overall sign of the pion exchange, for physical values of the masses, is much less 
model dependent than the previous chiral perturbation theory analysis seems to imply. In fact, if the 
7i-°^*<-y* form factor does not change its sign in the Euclidean range of integration over fcj, the overall 
sign is fixed even without knowledge of the form factor. This implies the same positive sign without 
use of the chiral limit, i.e. the same sign for exchanges of heavier pseudoscalars, J PC = where 
no large logarithms are present. Moreover, one can verify the same positive sign for exchanges by 
mesons with J = 1 ++ , 2 ^ with an additional assumption about dominance of one of the form 
factors. Exchanges with J PC = ++ , 1 h , 2 ++ give, however, contributions with a negative sign to 
ft HLbL unc [ cr similar assumptions, but they are much smaller. 

Next-to-leading terms in the large N c limit 

Now let us turn to the next -to-leading terms in 1/N C expansion. Generically these terms are due to 
two-particle exchanges in the HLbL amplitude, see the diagram in Fig. 2(b) with 7r + 7r~ substituted by 
any two meson states. What is specific about the charged pion loop is its strong chiral enhancement 
which is not just logarithmic but power-like in this case. In Eq. (jHJ) it is reflected in the term C2 m^/m 2 . 
The point-like pion loop calculation which gives a HLbL (7T7r) = —4.6 x 10~ 10 corresponds to C2 = 
—0.065. The rather small value of C2 can be contrasted with the one of the coefficient c\ which is not 
suppressed: c\ ~ 1.7. As we will see the smallness of ci is related to the fact that chiral perturbation 
theory does not work in this case. To see that this is indeed what happens is sufficient to compare 
the point like loop result with the model dependent calculations where form factors are introduced. 
Two known results, a HLbL (7T7r) = -(0.4 ±0.8) x I0" 10 [lHj and a HLbL (7T7r) = -(1.9 ±0.5) x I0~ 10 
[3 IB], show a 100% deviation from the point-like number. It means that the bulk of the contribution 
does not come from small virtual momenta ki and, therefore, chiral perturbation theory should not 
be applied. In other words, the term C3 in Eq. ([8]) with no chiral enhancement is comparable with 
C2("ip/w 2 .). It means that loops with heavier mesons should also be included. 
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Breaking of the chiral perturbation theory looks surprising at first sight. Indeed, the inverse chiral 
parameter rrvi/m 2 . rs 30 is much larger than N c = 3. What happens is that the leading terms in the 
chiral expansion are numerically suppressed, which makes chiral corrections governed not by m \ / m 2 
but rather by « 40 m 2 / m 2 . This can be checked analytically in the case of the HVP contribution 
to the muon anomaly. The charged pion loop is also enhanced in this case by a factor m 2 /m 2 but 
the relative chiral correction due to the pion electromagnetic radius (evaluated with a cutoff at m 2 
in the irir spectral function) is ~ 40 m 2 /m 2 ln(m p /2m 7r ). Of course, if the pion mass (together with 
the muon mass) would be, say, 5 times smaller than in our real world, the charged pion-loop would 
dominate both in the HVP and the HLbL contributions to the muon anomalous magnetic moment. 

In concluding this Section, we see that the l/N c expansion works reasonably well, so one can use 
one-particle exchanges for the HLbL amplitude. On the other hand, chiral enhancement factors are 
unreliable, so we cannot limit ourselves to the lightest Goldstone-like states, and this is the case both 
for the leading and next-to-leading order in the 1/N C expansion. 



3. Short Distance QCD Constraints. 

The most recent calculations of 

a HLbL in 

the literature [HIH1IE1IH] are all compatible with the QCD chiral 
constraints and large-iV c limit discussed above. They all incorporate the 7r°-exchange contribution 
modulated by 7r°7*7* form factors Tiki, k 2 ), correctly normalized to the ir° — ► 77 decay width. They 
differ, however, in the shape of the form factors, originating in different assumptions: vector meson 
dominance (VMD) in a specific form of Hidden Gauge Symmetry (HGS) in Refs. a different 

form of VMD in the extended Nambu-Jona-Lasinio model (ENJL) in Ref. [H] ; large-7V c models 
in Refs. [HH]; and on whether or not they satisfy the particular operator product expansion (OPE) 
constraint discussed in Ref. [5], upon which we next comment. 

Let us consider a specific kinematic configuration of the virtual photon momenta ki,k 2 , k% in the 
Euclidean domain. In the limit q = these momenta form a triangle, k± + k 2 + k% = 0, and we consider 
the configuration where one side of the triangle is much shorter than the others, k\ « k 2 ^> k 2 . 
When k\ ss k 2 3> m 2 we can apply the known operator product expansion for the product of two 
electromagnetic currents carrying hard moments k\ and k 2 , 

d 4 x 1 Jd 4 x 2 e- lk ^ k ^Mx 1 ) j p (x 2 ) = ^-e„ p5l k s Jd 4 ze-^- z f 5 (z) + o(j^) ■ (10) 

Here j§ = Y] q Q 2 <77 7 75<Z is the axial current where different flavors are weighted by squares of their 

electric charges and k = (fci — k 2 )/2 w k\ sa — k 2 . As illustrated in Fig.[3]this OPE reduces the HLbL 
amplitude, in the special kinematics under consideration, to the AVV triangle amplitude. 



q — 



YVs 





Figure 3: OPE relation between the HLbL scattering and the AVV triangle amplitude. 

There are a few things we can learn from the OPE relation in Eq. p0|) . The first one is that the 
pseudoscalar and pseudovector meson exchanges are dominant at large k\^- Indeed, only 0~ and 1 + 
states are coupled to the axial current. It also provides the asymptotic behavior of form factors at 
large k\ « k\. In particular, we see that the 7r°7*7* form factor T(k 2 ,k 2 ) goes as 1/k 2 and similar 
asymptotics hold for the axial-vector couplings. The relation in Eq. (|10[) does not imply that other 
mesons, like e.g. scalars, do not contribute to HLbL, it is just that their 7*7* form factors should fall 
off faster at large k\ 2 . 
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The AW triangle amplitude consists of two parts: the anomalous, longitudinal part and the 
non-anomalous, transverse one; we consider the chiral limit where m 2 — > 0. Because of the absence 
of both perturbative and non-perturbative corrections to the anomalous AVV triangle graph in the 
chiral limit, the pion pole description for the isovector part of the axial current works at all values 
of fcf connecting regions of soft and hard virtual momenta. This, in particular, implies the absence 
of a form factor F(0, fcf) in the vertex which contains the external magnetic field. At first sight, this 
conclusion seems somewhat puzzling because for non-vanishing external momentum q the form factor 
F{q 2 , fcf) certainly is attributed to the pion exchange. The answer is provided by the observation that 
this form factor enters not in the longitudinal anomalous part, but in the transverse part. It is for 
this reason that the axial anomaly is not corrected by the form factor. In the transverse part the form 
factor shows up together with the massless pion pole in the form 

^( g 2 ,fc 3 2 )-^(0,0) 

(h + q) 2 ■ [ ' 

At q = this combination contains no pion pole at fc 2 = . It means that the discussed piece conspires 
with the pseudovector exchange to produce the transverse result and in this sense becomes part of 
what could be called the pseudovector exchange. It provides the leading short-distance constraint 
for the pseudovector exchange. Contrary to the case of the longitudinal component, the transverse, 
non-anomalous part of the AVV triangle is, however, corrected non-perturbatively [10] lllj. 

Additional constraints on subleading terms in the .F(fc 2 ,fc 2 ) form factor, which were derived in 
Ref. [T2], are also taken into account in the calculation quoted in Ref. 

The large momentum behavior which singles out pseudoscalar and pseudovector exchanges is, how- 
ever, not sufficient to fix per se a unique model for the evaluation of a HLbL because the bulk of the 
integral in Eq. (J2|) comes from momenta ki of the order of an hadronic scale. However, the faster 
decreasing of exchanges other than pseudoscalar and pseudovector ones makes these contributions 
numerically smaller. Moreover, the importance of asymmetric momenta configurations with two mo- 
menta much larger than the third one was checked in [T3] numerically. This check is related to a 
question which we next discuss. 

There are other short-distance constraints than those associated with the particular kinematic 
configuration governed by the AVV triangle. At present, none of the light-by-light hadronic pa- 
rameterizations made so far in the literature can claim to satisfy fully all the QCD short-distance 
properties of the HLbL amplitude which is needed for the evaluation of Eq. Q. In fact, within the 
large-7V c framework, it has been shown [H] that, in general, for other than two-point functions and 
two-point functions with soft insertions, this requires the inclusion of an infinite number of narrow 
states. However, a numerical dominance of certain momenta configuration could help. In particular, 
in the model of Ref. [9] with a minimal set of pseudoscalar and pseudovector exchanges, the corrections 
due to additional constraints not satisfied in the model turn out to be quite small numerically. Note 
that in the frameworks of the ENJL model [5] the QCD short-distance constraints are accounted 
for by adding up the quark loop with virtual momenta larger than the cutoff scale of the model. 

4. Hadronic Model Calculations 

In the previous section we have mentioned a few models used for the calculations of a HLbL : HGS model 
in [H El [6] , ENJL model in [3 [8] , the pseudoscalar exchange only in pQ , the OPE based model of 
pseudoscalar and pseudovector exchanges in [9] . In order to compare different results it is convenient 
to separate the hadronic light -by-light contributions which are leading in the l/iV c -expansion from 
the non- leading ones [T5] . 

Contributions which are leading in the 1/N c expansion 

Among these contributions, the pseudoscalar meson exchanges which incorporate the 7r°, and to a 
lesser degree the rj and rj exchanges, are the dominant ones. As discussed above, there are good QCD 
theoretical reasons for that. In spite of the different definitions of the pseudoscalar meson exchanges 
and the associated choices of the T{k\, A: 2 ) form factors used in the various model calculations, there 
is a reasonable agreement among the final results, which we reproduce in Table 1. 



G 



Table 1: Contribution to a HLbL from 7r°, i] and rj' exchanges 



Result 


Reference 


(o c i 1 o\ y i n-10 

(8.3 ±0.6) x 1CT 10 
(8.3 ± 1.2) x 1CT 10 
(11.4 ± 1.0) x lO" 10 


[7l [81 

DO 



In fact, the agreement is better than this table shows. One should keep in mind that in the ENJL 
model (the first line) the momenta higher than a certain cutoff are accounted separately via quark 
loops while in the OPE based model these momenta are already included into the result (the last 
line in the Table 1). Assuming that the bulk of the quark loop contribution is associated with the 
pseudoscalar exchange channel one gets 10.7 x 10~ 10 in the ENJL model instead of 8.5 x 10 -10 . In the 
calculations quoted in the two other entries, the higher momenta were suppressed by an extra form 
factor in the soft photon vertex and no separate contribution was added to compensate for this. 

Closely related to pseudoscalar exchanges is the exchange by the pseudovectors. Both enter the 
axial-vector current implying relations between form factors (see the discussion of the triangle am- 
plitude in the previous section). Again, here the estimates in the literature differ by the shape of 
the form factors used for the ^7*7* and ^7*7 vertex. Different assumptions on hadronic mixing is 
another source of uncertainty. Although the contribution from axial-vector exchanges is found to be 
much smaller than the one from the Goldstonc-likc exchanges by all the authors, the central values, 
shown in Table 2, differ quite a lot. The authors of Ref. pj] attribute this to the influence of the OPE 

Table 2: Contribution to a HLbL from axial-vector exchanges 



Result 


Reference 


(0.25 ±0.10) x 10~ 10 
(0.17 ±0.10) x 10~ 10 
(2.2 ±0.5) x 10~ 10 


mum 
m 



constraint for the non-anomalous part of the AVV triangle amplitude, discussed above. Further study 
of the discrepancy in this channel is certainly needed. 

The scalar exchange contributions have only been taken into account in Refs. [TJ, [5] . In fact, within 
the framework of the ENJL model, these contributions are somewhat related to the constituent quark 
loop contribution. The result is: 

Contribution to a HLbL from Scalar exchanges [H] 
-(0.7 ±0.2) x 10~ 10 . 

It is much smaller than the contribution from the Goldstone-like exchanges and negative. In compar- 
ison with the pseudovector exchange, the magnitude for the scalar is a few times smaller than for the 
pseudovector in the OPE-based model but a few times larger in HGS and ENJL models. 

As we discussed in Section 2 there is some number of other C-even mesonic resonances in the 
mass interval 1-2 GeV, not accounted for in the ENJL model, which could contribute to a HLbL 
comparably to the contribution from scalars. These contributions are of both signs depending on 
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quantum numbers. At the moment we can only guess about their total effect. Thus, it seems reasonable 
to use the scalar exchange result rather as an estimate of error associated with these numerous 
contributions. 

Contributions which are subleading in the 1/N c expansion 

As we discussed in Section 2 the charge pion loop chirally enhanced as m} p jrn\ is a priori the dominant 
contribution in the subleading l/N c order. It occurs, however, that the chiral enhancement does not 
work and loops involving other heavier mesons can compete with the simple pion loop contribution. 

The dressed pion loop results are considerably smaller than the one for the point-like pion. They 
are presented in Table 3. The last line from Ref. [5] is not the result of a calculation. Strictly speaking 

Table 3: Contribution to a HLbL from a dressed pion loop 



Result 


Reference 


-(0.45 ±0.85) x 10- 10 
-(1.9 ±0.5) x 10~ 10 
(0± 1) x 10~ 10 


[HE] 
0! 




it represents an error estimate of the meson loop contributions subleading in l/A^-expansion. One 
can probably increase this error to cover the ENJL result in the second line. 

5. Numerical Conclusions and Outlook 

What final result can one give at present for the hadronic light-by-light contribution to the muon 
anomalous magnetic moment? It seems to us that, from the above considerations, it is fair to proceed 
as follows: 

Contribution to a HLbL from tt° , rj and rf exchanges 

Because of the effect of the OPE constraint discussed above, we suggest to take as central value the 
result of Ref. [9] with, however, the largest error quoted in Refs. Q3E]: 

a HLbL (7r , V , T)') = (11.4 ± 1.3) x 10- 10 . (12) 

Let us recall this central value is quite close to the one in the ENJL model when the short-distance 
quark loop contribution is added there. 

Contribution to <j HLbL from pseudovector exchanges 

The analysis made in Ref. [H] suggests that the errors in the first and second entries of Table 2 are 
likely to be underestimates. Raising their ±0.10 errors to ±1 puts the three numbers in agreement 
within one sigma. We suggest then as the best estimate at present 

a HLbL (pseudovectors) = (1.5 ± 1) x 10~ 10 . (13) 
Contribution to a HLbL from scalar exchanges 

The ENJL-model should give a good estimate for these contributions. We keep, therefore, the result of 
Ref. [TJ [5] with, however, a larger error which covers the effect of other unaccounted meson exchanges, 

a HLbL (scalars) = -(0.7 ± 0.7) x 10~ 10 . (14) 
Contribution to a HLbL from a dressed pion loop 
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Because of the instability of the results for the charged pion loop and unaccounted loops of other 
mesons, we suggest using the central value of the ENJL result but wit a larger error: 

a HLbL (7r-dressed loop) = -(1.9 ± 1.9) x 1(T 10 . (15) 

From these considerations, adding the errors in quadrature, as well as the small charm contribution 
in Eq. ([6]), we get 

<z HLbL = (10.5 ± 2.6) X 10" 10 , (16) 

as our final estimate. 

We wish to emphasize, however, that this is only what we consider to be our best estimate at 
present. In view of the proposed new — 2 experiment, it would be nice to have more independent 
calculations in order to make this estimate more robust. More experimental information on the decays 
7T° — * 77*, 7T° — * 7*7* and tt° — > e + e~ (with radiative corrections included) could also help to confirm 
the result of the main contribution in Eq. (fT^|) . 

More theoretical work is certainly needed for a better understanding of the other contributions 
which, although smaller than the one from pseudoscalar exchanges, have nevertheless large uncertain- 
ties. This refers, in particular, to pseudovector exchanges in Eq. (fT3|) but other C-even exchanges 
are also important. Experimental data on radiative decays and two-photon production of C-even 
resonances could be helpful. An evaluation of l/iV c -suppresed loop contributions present even a more 
difficult task. New approaches to the dressed pion loop contribution, in parallel with experimental 
information on the vertex 7r + 7r _ 7*7*, would be very welcome. Again, measurement of the two-photon 
processes like e + e~ — > e + e~ n + n~ could give some information on that vertex and help to reduce the 
model dependence and therefore the present uncertainty in Eq. (|15|) . 
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